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We evaluate numerically the momentum-space gluon propagator in the lattice Lan- 
dau gauge, for three- and four-dimensional pure 51/(2) lattice gauge theory. Al- 
though there are large finite-size effects, we always observe, in the limit of large 
lattice volumes, a gluon propagator decreasing in the infrared limit. This result 
can be interpreted in a straightforward way, by considering the proximity of the so- 
called first Gribov horizon in the infrared directions. We also consider the problem 
of discretization errors introduced by the lattice regularization, and their effect on 
the ultraviolet behavior of the gluon propagator. 

1 Introduction 

The gluon propagator is not an observable, since it is a gauge-dependent quan- 
tity. Nevertheless, the study of its infrared behavior provides us with a power- 
ful tool far increasing our understanding of confinement in non-Abelian gauge 
theories. Ill In fact, the infrared behavior of the gluon propagator can be di- 
rectly related to the behavior of the Wilson loop at large separations and to 
the existence of an area law.H 

The theoretical predictions for the behavior of the Landau-gauge gluon 
propagator in the infrared limit - obtained mostly by solving the gluon Dyson- 
Schwinger equation, approximately - range from a p""* singularity □ to a van- 
ishing propagator. □ A true nonperturbative investigation of these predictions 
is possible by performing Monte Carlo simulations of QCD on the lattice. 

2 The lattice setup 

Let us consider a standard Wilson action for SU (2) lattice gauge theory in d 
dimensions. The lattice gluon propagator (in momentum space) can be written 
as 
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where 

D,{k) ^ (3) 
Here V is the lattice volume, 

A^{k) EE A^{x) exp [27ri (fc • a; + (4) 

X 

and the lattice gluon field A^{x) is given by 

A,{x) ^ 1 [C/^(x) - {/t(x)] . (5) 

In order to fix the lattice Landau gauge we can look for a local minimum 
of the functional 

^u[9] ^^-^J2Y. [9(.^)U,ix)gHx + e,)] . (6) 

fl—l X 

In fact, if the configuration {[/^ (a;)} is a stationary point of the functional 
£u[9] then the lattice divergence of A^i(a;) is zero, i.e., 

d 

(V • A) (x) EE J2 iM^) " - e,)] = . (7) 

This is the lattice formulation of the usual (continuum) Landau gauge-fixing 
condition. Moreover, requiring this stationary point to be a minimum of the 
functional £u[g] implies that the transverse gauge-fixed configurations belong 
to the region delimited by the so-called first Gribov horizon, defined as 
the set of configurations for which the smallest non-trivial eigenvalue of the 
Faddeev-Popov operator is zero. 

Thus, when the lattice Landau gauge is imposed, theLphysical configuration 
space is restricted to the region and one can prove □ a rigorous inequality 
for the Fourier components of the gluon field From this inequality it 

follows that the region is bounded by a certain ellipsoid 8. This bound 
implies the proximity of the first Gribov horizon in infrared directions and 
the consequent suppression of the low-momentum components of the gauge 
field. This bound also causes a strong suppression of the gluon propagator in 
the infrared limit. In fact, Zwanziger proved □ that, in four dimensions and 
in the infinite- volume limit, the gluon propagator is less singular than p^^ in 
the infrared limit and that, very likely, it does vanish in this limit. A similar 



2 



result holds in three dimensions: one obtains that, in the infinite- volume limit, 
the gluon propagator must be less singular than p^^ as p —>■ and that, very 
likely, it vanishes in the infrared limit. We remark that these predictions for 
the gluon propagator are (3 -independent: in fact, they are derived only from 
the positiveness of the Faddeev-Popov operator when the lattice Landau gauge 
is imposed. 



Iff 



3 The infrared behavior of the gluon propagator 

We have studied the momentum-space gluon propagator D(k) in four0 and in 
threep dimensions. In both cases we have found that, if the lattice volume V is 
large enough, the gluon propagator is decreasing as the magnitude of the lattice 

[ d 1-^/^ 

momentum p{k) = 2 ^i^^ i'^^fj.) decreases, provided that p{k) is 

smaller than a value pdec- Also, the lattice volume at which this behavior 
for the gluon propagator starts to be observed increases with the coupling /3, 
i.e., finite-size effects are very large in the small-momenta sector. This makes 
practically unfeasible, with present computational resources, the numerical 
study of the infrared behavior of the gluon propagator in four dimensions and 
at large values of /3, and explains why a decreasing gluon propagator has b£^ 
observed only in the strong-coupling regime for the four-dimensional case.L 

In the three-dimensional case we get good scaling for values of f3 greater 
than 3.4, especially in the region of large momenta, where finite-size effects 
are negligible. We also see that the gluon propagator is decreasing for mo- 
menta p < Pdec, and that the value of pdec (in physical units) is practically 
/3-independent. From our data we obtain pd^c ~ 350 MeV. For the same set 
of data we also observe that the gluon propagator is less singular than p^^ 
in the infrared limit, in agreement with Zwanziger's prediction. We notice 
that the turnover momentum pto ~ 700 MeV is in good agreemfint with the 
result obtained recently in four dimensions for the SU{S) group. E3 Finally, the 
value D{0) of the gluon propagator at zero momentum decreases monotoni- 
cally as the lattice volume increases (see for example the case f3 — 5.0 in the 
three-dimensional case El). This suggests a finite value for D{0) in the infinite- 
volume limit, but it is not clear whether this value would be zero or a strictly 
positive constant. Therefore, the possibility of a zero value for D{0) in the 
infinite- volume limit is not ruled out. 



3.1 Conclusions 

The predictionBi of a gluon propagator decreasing for momenta p{k) < pdec is 
clearly verified numerically for several values of the coupling /3. In the three- 
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dimensional case these values range from the strong-coupling regime to the 
scaling region. Also, our data in the strong-coupling regime for the thrce- 
dimensionaljXase are in qualitative agreement with the results obtained in four 
dimensions. 13 This strongly suggests to us that a similar analogy will hold — 
in the limit of large lattice volumes — for couplings (3 in the scaling region, 
leading to an infrared-suppressed gluon propagator also in the four-dimensional 
case. 

4 Discretization effects 

The definition of the lattice gluon field given in Eq. (||) is only one of the 
possible lattice discretizations of Afj,{x), i.e., we can consider several possible 
definitionsT,leading to discretization errors of different orders. For example, we 
can write 




(8) 



2i 




(9) 



(^) 



(10) 



8i 



If we set 




we obtain that A^i}\x), A'j^\x) and A'j^\x) are equal to agQ a-A{x) plus terms 
of order gg- We can also consider 



4'^)(x)^ [4 4)(:r) - 42)(x)J /3 (12) 
45)(x)^[l6 4)(x)-43)(x)] /15 (13) 
4'^)(x)^ 4 42)(x) - Al'\x)] /3 (14) 



and 



A<^J\x) ^ 64 4) (x) - 204^) (a:) + A<^^\x) /45. (15) 
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Table 1: The parameter Z and the jA.o.i. for the best fit to the ultraviolet behavior 
Zqq/ |^4p^{fc)j. We consider different lattice discretizations of the gluon propagator, and 

indicate with Dy|(A:) the D^^'>{k) propagator evaluated using tadpole-improved links Ui_i(x). 



Type of 
propagator 


I3-- 


= 2.7 




= 10 


Z 




Z 


xVd.o.f 




1.2965 
1.6052 
0.9141 
1.4012 


1.36 
1.36 
1.25 
1.25 


1.0261 
1.0904 
0.9295 
1.0496 


0.121 
0.121 
0.133 
0.133 



The definitions A^^^ and A^^^ were recently considered by Giusti et al. E3 
They find that the corresponding gluon propagators are equal modulo a con- 
stant factor. We have performed a similar studyEll at several values of (3 and for 
lattice volumes V — 8'^ and 12'^, evaluating D^*^ {k) using the different defini- 
tions of the gluon field yllt''' given above. We obtain, in all cases, that the seven 
propagators D^'^\k) are equal modulo a constant factor. We notice that this 
proportionality constant between different discretizations of the gluon propaga- 
tor may be explained as a tadpole renormalization. In fact, let us consider the 
tadpole-improved link U^{x) = Uf^{x) / uq, where uq is the mean link in Landau 
gauge. Then D'^^\k) gets multiplied by a factor Uq^, and D^'^\k) by Uq^- For 
example, for V = 12'* and f3 = 2.2 we have Uq = 0.68428(8). In this case, using 
tadpole-improved operators, the discrepancy D''^\k) / D^'^^ [k) is reduced from 
1.833(4) to 1.254(3). Similarly, at 13 = 2.3 [and 2.7], we obtain ul = 0.7222(1) 
[respectively Uq = 0.8077(3)] and the discrepancy D'^^'' {k) / D'^'^\k) is reduced 
from 1.701(4) [respectively 1.425(3)] to 1.228(3) [respectively 1.151(2)]. As 
expected, the discrepancy between different discretizations decreases as (3 in- 
creases. 

It is also interesting to compare the data obtained for the gluon propa- 
gator at large momenta with the ultraviolet behavior g^/ [4p^(fc)] , predicted 
by perturbation theory at zeroth order. To this end, we consider different 
discretizations of the gluon propagator and we fit the data corresponding to 
p^(fc) > 3 using the function Zg^/ [4p^(A;)]. In Tab. |l] we report the results 
obtained for the lattice volume V = 12** at /3 = 2.7 and /3 = 10. Again, the 
discrepancy between D'^^\k) and Z?(^)(fc) is reduced when tadpole-improved 
operators D!^\{k) are used. Also, as expected, the value of Z gets closer to 1 
as /3 increases. 
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Finally, let us notice that, in finite-temperature QCD, the long-distance 
behavior of the gluon propagator is directly related to the electric and magnetic 
screening lengths, and that these screening masses are invariant under rescaling 
of the propagators by a constant factor. In particular, D[^\k) + D2 \k) 
[resaectively D^'^^k)] is related to the gluon propagator used by Karsch et 
aLll3 for the evaluation of the magnetic [respectively electric] screening mass. 
It has been checkedliil that, also at finite temperature, different discretizations 
for the gluon propagator are equal modulo a constant value, i.e., the screening 
masses are independent of the discretization 13^*' (fc). 
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